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Abstract 

^ Using the general expressions for level shifts obtained from the master equation for a small 

I system interacting with a large one considered as a reservoir, we calculate the dispersive po- 

tentials between an atom and a wall in the dipole approximation. We analyze in detail the 
particular case of a two-level atom in the presence of a perfectly conducting wall. We study 
. the van der Waals as well as the resonant interactions. All distance regimes as well as the high 

■ and low temperature regimes are considered. We show that the Casimir-Polder interaction can 

not be considered as a direct result of the vacuum fluctuations only. Concerning the interac- 
tion between the atom and the wall at high temperature, which show that a saturation of the 
potential for all distances occurs. This saturated potential coincides exactly with that obtained 
in the London- van der Waals limit. 
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1 Introduction 



The subject of dispersive interactions is an old problem related to the intermolecular forces whose 
discussion first eppeared in the van der Waals state equation for real gases in 1873 pQ. In this 
equation there is a term proportional to 1/v 2 , where v is the specific volume of the gas, responsible 
q-i for a long-range interaction between the molecules. This interaction provides a mechanism for the 

phase transitions between the gas and liquid states (and vice-versa) of the substance. However, only 
in the beginning of the 20th century, these interactions were understood [5] ■ 

Approximately a hundred years ago, M. Reinganum [3] suggested that van der Waals forces 
could be explained by the electromagnetic interaction between the electric dipoles associated to the 
molecules. In the twenties, mainly due to W.H. Keesom and P. Debye works [H[S], Reinganum's 
insight showed to be very fruitful, making possible the explanation of thevan der Waals interaction 
between two polar molecules (orientation forces) and between a polar molecule and a nonpolar one 
(induced forces). 

However, the explanation for the force between neutral and non-polar molecules, like noble gases 
for example, had to wait for the advent of Quantum Mechanics. In 1930, F. London [5] obtained 
the interaction potential between two neutral hydrogen atoms in their respective ground states and 
separated by a distance much larger than Bohr radius in the static electric dipole limit, Ao S> r 3> ao, 
where Ao = 2irc/ujQ is the wavelength of the main atomic transition relevant to the interaction and 
ao is the Bohr radius. London's result is V(r) ~ — 6.5e 2 ag/r 6 ~ — 3Uuj al / ir 6 . 

As one can see, this interaction is a genuine quantum effect (proportional to H) and depends only 
on the atomic polarizabilities. Since the polarizability is intimately connected to dispersive part of 
the electric permittivity of the substance, this kind of interaction is called dispersive van der Waals 
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interaction [7]. Then, dispersive van der Waals forces are long-range electromagnetic forces between 
non-polar and neutral molecules in their ground states. 

In 1948, Casimir and Polder [Sj considered for the first time the influence of retardation effects 
on the van der Waals forces between two atoms as well as on the force between an atom and a 
perfectly conducting wall. They showed that in the retarded limit the interaction falls as 1/r 7 for 
two atoms (in contrast to London's 1/r 6 result) and as 1/r 4 for an atom and a perfectly conducting 
wall, in contrast with the short distance limit (that falls as 1/r 3 ) [9]. 

In 1956, E.M. Lifshitz developed a general theory of dispersive van der Waals forces between 
dielectric macroscopic bodies using a non-perturbative approach [10j . because the perturbative ex- 
pansion for many-body interaction breaks down due to the non-additivity of the van der Waals 
forces. Lifshitz derived a powerful expression for the force at finite temperature between two semi- 
infinite dispersive media characterized by an electric permittivity and separated by a slab of any 
other dispersive medium. He was able to derive and predict several results, like the variation of the 
thickness of thin superfluid helium films in a remarkable agreement with experiments |11] . He also 
showed that the Casimir-Polder force is a limiting case when one of the media is sufficiently dilute 
such that the force between the slabs may be obtained by direct pair-wise integration of a single 
atom- wall interaction [T2] . 

Since then a wide knowledge about the behavior amd nature of dispersive forces has been 
achieved. In the following decades, many works have been done about this subject, like those 
done by A.D. McLachlan [T3], in the 60's, where the recently compiled linear response theory devel- 
oped by R. Kubo [14] was applied. In 1970, G. Feinberg and J. Sucher [15] treated the electric and 
magnetic contributions to van der Waals interaction in the same foot. In the last 60 's decade and 
during the 70's it was analyzed the interaction between excited atom and a wall [T5], [T7] [T5] . The 
first of these works, made by H. Morawitz, considered the level and frequency shifts of an excited 
atom in the presence of a perfectly conducting wall using the image method, which led him to the 
discovery of the resonant dispersive interaction, that we will discuss later. In the 80's and 90's the 
non-additivity character of the van der Waals forces was exhaustively exploited pjH [20] ; very good 
works concerning level shifts of atoms and dispersive interactions in cavities [21] (22] [23l [24] [25] and 
the thermal contributions to the Casimir-Polder interaction between atoms and finite systems (like 
cylinders or spheres), in a microscopic approach, was studied [2"6l |2"T1 12"8"] . 

More recently the interest for van der Waals and Casimir-Polder forces have been raised mainly 
due to better experimental techniques and the development of nanosciences. The influence of these 
forces in the stability of Bose-Einstein condensates [29] [30] and on factoring and doping carbon 
nanotubes are branches of great activity nowadays [311 E21 [33] . 

In this paper we make the sequence of the work began in a previous one [34] . where we studied 
the van der Waals interaction between a two-level atom and a perfectly conducting wall using the 
density matrix formalism [35] , Here, we come back to this problem using the same formalism, 
but studying also the excited state contributions to the interaction, both in vacuum and thermal 
states of the electromagnetic field. We give simple interpretations to the ground state contribution 
(van der Waals interaction) and the excited state contribution (resonant interaction) in terms of 
the non-resonant and resonant parts of the atomic polarizability. Considering the thermal state 
of the electromagnetic field, we analyze both low and hight temperature limits, Two® <C fcgT and 
Hujq 3> ksT respectively, and show that in the last limit the interaction between the two-level system 
and the perfectly conducting wall is temperature independent and behaves exactly like that in the 
London-van der Waals limit for vacuum contribution. 

2 Level shifts and exchange energy rates 

In this section we obtain the level shifts and the exchange energy rates of a multi-level system 
interacting with the electromagnetic field in the dipole approximation. In the first case, if we consider 
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the boundary conditions imposed to the field by the presence of a given wall, we can extract the 
position dependent part, which gives, by derivation, the atom-wall dispersive forces. 



2.1 General expressions 

Our starting point is the master equation for a small system S (the atom) which interacts weakly 
with a large one TZ, that may be considered as a reservoir (quantized electromagnetic field), by the 
interaction hamiltonian 

V = -J2s j R j , (1) 

3 

where Sj and Rj are observables associated to S and TZ respectively. In the Schrodinger picture, 
the master equation may be written as [351 1361 137j 



^Pab (*) = -^abPab (*) + E helped (*) » 



(2) 
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where pf fc = (a\ps\b) is the matrix element of the density operator p$, associated to the system 5, 
between the energy eigenstates \a) and \b) with eigenvalues E a and Eb of the unperturbed hamil- 
tonian operator Hs of S (Hs\a) — E a \a))\ S° ab = (a\Sj\b) is the matrix element of the observable 
Sj and uj a b = {E a — E^) /Ti is the transition frequency between the states \a) and \b). The function 
9jk ( r ) that appears in equation Q is defined as 



■9,-fc 



to = [flg (-r)] * = Tr fl (r) (0)] = £p„ £ R^R^e^ , 



(4) 



where = (fj,\Rj\v) is the matrix element of the observable Rj between the energy eigenstates 
\p) and \v) with energy eigenvalues E^ and E v of the unperturbed hamiltonian operator Hr of TZ 
(Hji\/i) = E^l/j,)); uj^ v = (E^ — Ey) /K is the transition frequency between the states and \v) 
and pr is the density matrix associated to TZ, considered constant in time and diagonal in the {\p}} 
base-ket, so that 



pr. = ^2p^\p)(p\ > 



(5) 



where p M is the statistical weight of the state \p) for a given ensemble. 

Making a brief resume, equations were obtained from the equation of the time evolution of 
the total density matrix p (t) ~ ps (t) <%> p^, 
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[p(t),H] 



where H = Hs + Hr + V is the hamiltonian of the total system S + 1Z. Then, one has made 
implicitly the assumption that there are two very different time scales in the evolution of the total 
system S + 1Z: the characteristic time Tg in which the average values of the observables of S change 



significantly and the characteristic time t c , which is, crudely speaking, the width of g?L (r) and 



measures the time of the fluctuations of the reservoir observables. 
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The approximations used in derivation of the equations refered above are based on the following 
conditions, 

r c « At « T s , (6) 

where At is the time interval that enters in the calculation of the time derivative in equation (|2j). 
Last condition tells us that the master equation is a coarse grained rate, since rapid variations of ps 
that occur in times of order of t c are smoothed in the interval At; since we are interest in times of 
the order of Tg, this smoothing remains a good aproximation once condition ([6]) is satisfied. Another 
suplementary condition for obtaining equation ((2|), which makes the coefficients T a \, c d independent 
of At, is that frequency difference between the elements of the density matrix should be very small 
compared to 1/At, so that 

l^ac - ^bd|At < 1 . 

Last condition is called secular condition and the sum in equation ^ is only over secular terms 
(terms that satisfy the secular condition). 

Using two new assumptions, one may conveniently split equation @ into two other ones: an 
equation for the diagonal elements pf a (populations) and another for the non-diagonal elements pf b , 
with a b (coherences). For the populations, the assumption is that there is no coherence p^ d in 
the sum in equation ([2]) with too low frequency, so that u c dTs S> 1. This implies that all coherences 
are not secular terms and the master equation for the populations pf a may be put into the form 

f t pL(t) = r c ^ a -p? a (i)r ^ c ) , (7) 

c 

r c ^a = ^E^E^^M! 2 ^^-^-^) ' ( 8 ) 

where V is the interaction hamiltonian given by |T]). The quantity r c ^ a may be interpreted as the 
transition rate probability between the states \c) and \a) as a result of the interaction of the system 
with the reservoir. 

In order to obtain the equation for coherences we make the second assumption, namely, we 
consider only the non-degenerate case: the frequency u; ob differ from all other frequencies uj ci i by a 
quantity of order (or larger than) 1/At, that is, \uj ac — ujbd\ ~ 1/At 3> I/T5. Then, the only secular 
term in equation @ is that which couples the coherence p^ b with itself, so that the equation for this 
coherence may be written as 

^Pab (*) = -* Kh + A Qb ) p s ab (t) - r QbP f b (t) , a ? b (9) 

where 



A ab = A a -A b (10) 

r ab = L^ + rr d - + rr d ' (11) 

^ = AE^EE^ l(M ' n|m ^ )|2 a*) 

fj. V J V J 

Kf = -^E^E^^I^M^M^) (13) 

rr d - = pE^EEi^"i^^)i 2(5 (^-^o = ^E r ™^ ( 14 ) 

H v j j 



and the symbol V in equation (|12p represents the Cauchy principal value. 
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The integration of the equation © is a straightforward task and gives 



pf h W=pf b (0)e- r -*e-^- +A -) t . 



from which one may conclude that the coherence p^ b , as a result of its interaction with the reservoir, 
oscillate with a frequency shifted by A a ;, from its free value Lo a b and decays exponentially with a 
characteristic time given by l/T a b. From equations (JTHHJ) and (fT2|) one may obtain the exchange 
energy rates between the system and the reservoir and the energy level shifts of the system, respec- 
tively, since frequency shifts result from levels shifts. Recall that SE a = hA a is the shift in the 
energy of the state \a) (and equivalently for \b)). 

For energy rates, considering the time derivative of the average value of the hamiltonian of the 
system in Schrodinger and Heisenberg pictures, one will find 
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dt 



— ^ {Eb — E a ) r a _,6 = ^ TlUJbaFa- 



(15) 



which may be interpreted as a net rate of changing the average value of the energy of the system S 
when it is in the state \a) or, in other words, the exchange energy rate between S and 1Z for the state 
| a). It is possible to put last result into a more intuitive form that clearly exhibits the roles played 
by S and TZ. Returning to equation (j4|) , one may show that real and imaginary parts of gj| (r) are 
related to symmetric correlation function (r) and the linear susceptibility xf k ( r ) E3 EE] 



Cf k (r) = Re [gf k (r)] = ^ p M ^ R^R^ cos (w^r) , 



k (r) = 76 (r) Im [gf k (r)] = - - ]T p, £ R^R^Q (r) sin (uj^t) 



(16) 
(17) 



where (r) is the step function. In the frequency space, we have 
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(18) 

(19) 
(20) 
(21) 



The symmetric correlation function given by ()18j) may be interpreted as a dispersion of the observ- 
ables Rj and R k arround a frequency w. It is a measure of the fluctuations of the reservoir dynamical 
variables. The real part of the susceptibility defined in equation (|20[) is called the dispersive or re- 
active part and is related to the polarization of a system by an external perturbation. On the other 
hand, the imaginary part is directly responsible for the absorption and dissipation of energy by the 
system and, for this reason is called dissipative part of susceptibility [39] - 

Now, let us come back to the equations (fT2|) and (fT5|) . Using the quantities defined in equations 
(|18H21[) . the energy level shifts and the exchange energy rates may be cast into the form 
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where 
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(24) 
(25) 
(26) 
(27) 

(28) 

(29) 
(30) 



are the symmetric correlation function and the dispersive and dissipative part of the susceptibility 
of the system 5 in the state |o) respectively. Equations (|24II27|) are the general expressions for the 
level shifts and the exchange energy rates of a small system interacting with a reservoir. 

The physical interpretation of equations (|24H25|) is simple. The former gives the contribution 
to the level shift due to polarization of the system S by fluctuations of the reservoir (fr) and the 
latter gives the contribution of the polarization of the reservoir by fluctuations of the system, or 
the reservoir reaction (rr) contribution. This last case is analogous to the radiation reaction of an 
accelerated charge in an external field. 

A similar interpretation may be done for equations (|26H27p . The first of them may be understood 
as the power absorbed by the system S from the reservoir fluctuations and the second represents the 
power dissipated by the system or, equivalently, the power lost to the reservoir. 



2.2 Dipole interacting with the radiation field 

Let us consider a small neutral but polarizable system which interacts weakly with the radiation 
field. Assuming the dipole approximation and considering the total quantized electric field operator 
at the positon x of the center of mass of the system as a sum over all possible modes, the interacting 
hamiltonian takes the form 

V = -d • E (x, t) = eJ2 £ (sj7£a W x + h.c.) , (31) 

kA j 

where d = — er is the dipole moment operator of the system, — e is the electrical charge, r = 
(xx,X2,Xs) is the position operator of the charge and the functions fkA (x) = ijf^ x (x) (we are using 
Einstein's convention here), where Xj is the unitary vector in the direction of component Xj of r, 
carry the information about the boundary conditions and possible source contributions. The index 
j runs over cartesian components (j = 1, 2, 3) and the index kA specifies a field mode, where k is its 
wave- vector and A is the corresponding polarization. Operators OkA and aj^ are the annihilation and 
creation operators, respectively, of one photon in the mode kA and satisfies the usual commutation 
relations 

[flkAjOk'A'] = [ a ikA' a k'A'] = J [ a kA,a k ,,v] = ^kk'^AA' ■ (32) 
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Choosing the coordinate axis in a manner such that the susceptibility and the correlation function 
of the system are diagonal we have 
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and 



2e 2 
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(33) 
(34) 
(35) 

(36) 
(37) 
(38) 

(39) 



is the static polarizability of the system in the direction Xj between the states \a) and \b). 

Using last equations in (|24H27[) , making = \rik\) in equation (which means a Fock state 
with n photons in the mode kA) and performing the integration on u>, we obtain, after some calcu- 
lations, 
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(40) 
(41) 
(42) 
(43) 

(44) 

(45) 



where k = uju/c, k ab = u a b/c and (nkA) is the average number of photons in the mode kA. 

Equations (|40H43[) are nothing more that equations (|24H27[) applied to a neutral but polarizable 
system interacting with the electromagnetic held. The validity of these equations, once assumed 
the dipole approximation, is determined by the validity of condition ([6]). This implies that the field 
correlation function in time, 



C l W = 2 E 1/kA (*) I 2 (W) + I) cos^r) , 



kA 
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obtained from equation (|16[) . must have all its characteristic frequencies much larger than any 
characteristic frequency of the system, so that the correlation time r c obeys the condition: t c <C 1/uiq, 
where ujq is the largest characteristic frequency of the system. 

Equations (|4H)l and (|4"Tj) . give the (rr) contributions to level shifts and energy rates of the system 
and equations (|42[) and (|43|) give the {fr) contributions. We see that (rr) contributions do not depend 
on (rikA)) which carries all information about the field state. This information is concentrated on 
{fr) contributions. Of course, for (n^x) = 0, there is still a residual term which can be associated 
with the contribution of the vacuum fluctuations of the field. 

Functions a'J^ (k) and cx^j^ {k) appearing in the {fr) terms are nothing more that the dis- 
persive and dissipative parts of the susceptibility of the system. In the (rr) contributions, however, 
we have the functions a'Jj "* (fc) and a"j ^ (fc), which do not have a direct interpretation like their 
corresponding "plus" functions. The difference is only in the sign of the k — kb a terms for the first 
and k + kb a for the second. But this fact will be crucial to determine the behavior of the dispersive 
interactions: it is precisely this difference of sign which makes possible the existence of the dispersive 
van der Waals interacion. This fact will become clear when we discuss the two-level system. 



3 Calculation of the dispersive forces 



In this section we shall apply the expressions deduced in the previous section for level shifts and 
energy rates of a polarizable system interacting with the electromagnetic field. We shall discuss a 
two-level atom, considered as an isotropic system, in the presence of a perfectly conducting wall. 

Let us consider two parallel square plates with sides L fixed at z = and at z = £, where L>(, 
and a small polarizable system on the OZ axes at position z that satisfies the condition < z -C I. 

Hence, for a given wave-vector k, the quantized electric field (without sources) in the region 
< z < £ is given by [3D]: 
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with n being a non-negative integer number. The mode with n = has an additional factor l/y/2 
not shown. Using this expression for the field in the limit £ — * oo, we shall be able to calculate the 
desired dispersive potentials between the atom and the single wall at z = 0. 



3.1 Two-level system 

First of all, let us assume isotropy, such that the polarizability given by equation (|39[) is independent 
of the space direction, that is, 

2e 2 

a le = a ge = -Ueg = 17=— | (ff|r|e) | 2 = a , (51) 
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where \g) and |e) are the ground and excited states of the two-level system under consideration with 
energies E g and E e , and ujq is the transition frequency, given by loq = fcoc = (E e — E g ) /Ti. 
Hence, equations (|4"0"f and (j4"2")) for level shifts take the form 
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(52) 
(53) 
(54) 



Note that, while (rr) contribution for level shifts of \g) and |e) states are exactly the same, (fr) 
contribution have the same magnitude but opposite signs. 

Using the expression (|46p in equations (|52M53p . taking the limit £ — > oo and considering a con- 
tinuous spectra for the field, we are faced with integrals of the form 
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where A > is a real parameter. For an analytical function / satisfying the condition 
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in all complex plane, it is possible to write 
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Hence, with the aid of last result, the z-dependent part of equations (|52j|53p , which gives the dis- 
persive potentials between the atom and the wall, when the average number of photons per mode is 
independ of the polarization, (n^x) = (rik), takes the form 
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(57) 
(58) 
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(60) 

(61) 

(62) 

(63) 
(64) 
(65) 

(66) 
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where (z) and Vj/ r (z) are the (rr) and {fr) terms of the vacuum contribution to the interaction 
and 7 is the Euler-Mascheronni constant. The term V/S {z) that carries the information about field 
state (thermal state, for example) and does not contribute when there is no photons in any mode of 
the field. 

Exploring this last case, that is, when (rik) — 0, which means that the field is in its vacuum 
state, we have 

V g {z,0) = V rr (z) + V Q fr (z) = V (z) , (67) 
V e {z,0) = V rr {z)-V fr {z) = 2V rr {z)-V {z), (68) 

for the ground state and excited state vacuum contributions to dispersive potentials, respectively. 
Analysing the small and large distance limits, one may show that 



For koz <C 1 (short distance limit): 



V e (z,0)^V g (z,0) = -^^ + O(z-i) . (69) 



For koz 3> 1 (large distance limit): 



V g {z,0) = -^+0{z- e ) ) (70) 

She ao + , 4 cos {2koz) 

7 + ncankr. 

8tt z 4 2k z 



V e {z,0) ~ — — + hca k — ^— — . (71) 



Equation (|69[) can be recognized as the interaction between the atom and the wall in the London- 
van der Waals limit [3]. In this limit, both excited and ground state potentials coincides and are 
mainly given by the (rr) contribution (only radiation reaction is important) as one can easily verify 
from equation (|59[) when k$z <C 1. 

The equation ([70]) . which gives the ground state potential at large distance, is the well known 
Casimir-Polder potential [8]. However, by equation (f7Tj) . we see that in this limit the excited 
state potential does not coincides with the ground state potential anymore. Apart from a Casimir- 
Polder term with opposite sign, there is an oscillating one that falls as 1/z [HI HPJ and since the 
Casimir-Polder term falls as 1/z 4 , the oscillatory term dominates over all large distance limit and 
the interaction presents an infite number of potential wells. 

In order to get a better understanding about last results, let us come back to equations (|52ti53|) . 
Making (n^x) = 0, the ground and excited state level shifts of the atom take the form 



SE, 



2 ^ k + k 
kA ' 



SE e = 5E? + SET = |fkA W |2 ■ (73) 

2 ^ k-kn 
kA u 

From last equations we clearly see that the potential associated to the ground state depends only 
on the non-resonant part of the atomic polarizability, while the potential associated to the excited 
state depends only on the resonant part. Since van der Waals forces are due to level shifts of the 
ground states of the interacting atoms [41) it may be called dispersive non-resonant interaction. On 
the other hand, for obvious reasons, the excited state potential may be called as dispersive resonant 
interaction. 

Last equations permit us also to understand why resonant interaction is very much stronger than 
van der Waals interaction at large distances and practically equal to this last at short distances. 
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The position dependent part of the modulus square of the amplitute of a field mode (subject to 
boundary conditions) which appears in the sums of the equations (|72M73[) . is usualy an oscillating 
function of products of characteristic distances and the frequency associated to the mode. In our 
case the characteristic distance is z, which means that for high frequencies such that k > 2tt/z, the 
oscillatory behavior of the summands leads to a small contribution to the total sum compared to 
the contribution comming from the low frequencies (k < 2tt/z). 

Taking into account that the number of modes per frequency is proportional to k 2 , at short 
distances virtual photons with high frequencies (A; ~ 2tt/z 3> ko) are much more important to the 
interaction than those with low frequencies (k <C ko), so that ko may be neglected in denominators 
of the sums in equations (|72H73p and the level shifts will be approximately the same. However, 
as the distance between the atom and wall becomes large, smaller and smaller frequencies become 
important to the interaction. For the ground state contribution, this leads to a smooth variation 
of the interaction. However, for the excited state contribution, when the distance becomes of the 
order 2~k jko, there is a strong increasing of the interaction, so the summand becomes singular at the 
frequency loq. At very large distances, small frequencies become more important and both ground 
and excited state potentials become weak, but the excited state potencial still remains very much 
stronger than the ground state potential. This happens because though the oscillations in summands 
tend to cancel contributions of frequencies 2it/z < k <C ko, the density number of modes is too small 
at very low frequencies, which still makes the principal value in (|T3[) the main contribution. 

One can get a simple interpretation of this fact. Looking at equation (fTTj) more carefully we 
easily see that a maximum and a minimum of the potential occur, respectively, at 



where Ao = 2nc/u)o is the transition wavelength and n > 1 is an integer (of course, these approxi- 
mations will be better for higher values of n) . Equations (|74|) are nothing more than the conditions 
for the existance of stationary waves with wavelength Ao in the "cavity" defined by the atom and 
the wall: the former is the resonant condition for a one-dimensional cavity with the "closed ends" 
and the latter is the resonant condition for a cavity with one "closed end" and one "opened end" . 
Then, it is comprehensible why the interacting potentials have the maximum magnitude or, in other 
words, a maximum response at that distances and not at others. 

In Figure Q] we show the exact potentials given by equations (|67H68j) and the (rr) and (fr) 
contributions given by equations (|59H60p . As we pointed out previously, the ground and excited 
state potentials are approximately the same as the London-van tier Waals potential (69) at short 
distances. At large distances, the excited state potencial has an oscillatory behavior while the ground 
state potential is given by the Casimir-Polder potential ([70]) . 

We saw that at short distances the radiation-reaction on the polarizable system is the dominant 
contribution to the interaction, a well established result in literature [16j ED [22] ■ However, Casimir- 
Polder potential is usually interpreted as a stark shift due to the field fluctuations mainly [2"2] . 
It is a common thing the use of expressions similar to equation (|53[) . for (n^x) — 0, in order to 
calculate the Casimir-Polder potential only replacing the atomic polarizability a' + (k) by its static 
value a' + (0) = ao in the summation [42]. This procedure can be justified if we invoke the argument 
that frequencies smaller than k ~ 2tt/ z <C fco give the main contribution to the interaction, so 
that the substitution of a' + (k) by ao causes a little error but gives the correct leading term of the 
potential. 

Now, despite the success of this procedure, this seems to hide the true physical content behind it. 
As one can see from Figure[T] both (rr) and (fr) contributions are of the same order of magnitude for 
distances z > 0.2Ao, but have opposite signs. For this reason, the excited state interaction given by 
the equation (|68|) is about two times (rr) contribution in the above mentioned interval. However, for 
the ground state potential given by the equation (|67[) , the (rr) and (fr) contributions almost cancel 



n 



— (maximum) , 




(minimum) , 



(74) 
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Figure 1: Dispersive potentials associated to the ground (van der Waals interaction) and excited 
(resonant interaction) states with corresponding (fr) and (rr) contributions. The total spontaneous 
emission rate as function of z, T spt (z), is also shown. The potentials are in units of hcaok^ and the 
spontaneous rate is in units of 2caokg, which is the value for the spontaneous rate in the free space. 



to each other, leading to a small potential value which is precisely the van der Waals interaction, 
which reduces to Casimir- Polder potential for z > Xo. Then, both vacuum field fluctuations and 
radiation-reaction play a crucial role for the interaction at large distances and the interpretation of 
the Casimir-Polder interaction as a shift due only to field fluctuations is not complete. 

In order to conclude this section, let us make a comment about the exchange energy rates for 
the two-level system. 

From equations (j¥T|) and ([4"3"|) , one may easily show that Q g = and 

Q e = -7ra fc 2 c^|f kA (x)| 2 (5(fc-fc ) , (75) 

kA 

where we considered the field in its vacuum state. For the geometry treated here, namely, with 
two-level atom at a distance z from an infinite perfectly conducting wall, last equation leads to [42] 

Q e = -hcu T spt (z) , r spt (z) = T s e tg [l-G (2k z) ] , (76) 

where r sp t (z) is the total spontaneous emission rate of the atom in the presence of the wall and 

Tlt g ^2ca k^—\(e\r\g)\W 

is the spontaneous emission for the atom in the empty space, being a — e 2 /hc ~ 1/137 the fine 
structure constant. Then, the life time of the excited state of the atom is of the order l/r spt (z), so 
that the resonant interaction is unstable and at a finite time (generally about 10 _7 s in the visible 
region [13]) the van der Waals interaction takes place. 

Equation (|75j) for r spt (z) is also plotted in the Figure [TJ The behavior of the spontaneous 
emission is similar to the behavior of the resonant potential. However, we should note that there 
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is a difference of phase between the potential and the spontaneous emission of approximately tt/2. 
Correspondimg to a maximum or a minimum of the potential, the spontaneous emission is exactly 
the same as if the atom was in the empty space. Somehow this could have been anticipated once for 
these cases, from conditions (ITU) , there is a node on the wall. 

A final comment is in order here: the fact that the potential is in quadrature with the spontaneous 
cmisson, may be understood if we remember that these quantities are related to real and imaginary 
parts, respectively, of the susceptibilities (of the atom and the field) and these quantities are, of 
course, in quadrature. 



3.2 Thermal corrections to the interaction 

Let us now turn our attention to the thermal corrections to dispersive interactions by considering 
both ground and excited states of the atom. There are many works in the literature treating thermal 
corrections to dispersive van der Waals forces, like [TUl [13 [23 to mention just a few. However, we 
have not found works that have analyzed also thermal corrections to the resonant interaction or 
the high temperature limit for dispersive interaction between the two level system and the perfectly 
conducting wall. This is our aim in this section. 

From equations (|57ti58[l . one may see that the only term that accounts for corrections to the 
interaction due to non-trivial states of the field is V/T (z), which is given by equation (|61|) . For a 
thermal state of the field at temperature T, we write 

v$> W = * W - v f G {2kz) dk ' (77) 

where At = hc/ksT is the thermal length, which defines the length scale beyond which thermal 
contributions pass to dominate over vacuum contribution of the van der Waals interaction. At room 
temperature, At — 7.63/im. 

In the low temperature limit, IcqXt S> 1, we have already caculated the last integral in both 
small and large distance limits in a previous work |34j , which for thermal corrections are defined by 
conditions z <C At and z > At respectively. 

For z <C At, we have 

VH^(T)-M^V, (78) 

where C(T) — 27r 3 ?Lca:o/45A T is independent of z and does not contribute to the force. The above 
result, that represents a small correction for the vacuum contribution, fits well the exact potential 
given by ([77| within an error smaller than 3% for z < 0.10 At and, since it does not depend explicitly 
on Aq, it is valid for both London- van der Waals and Casimir-Polder limitfl 
For z ~ At or larger, we have 

V (z) + V T (z) — VLif (z, T) = — —j— , (79) 

4 z 6 

which is exactly the Lifshitz asymptotic result [10j . which we hereafter call only Lifshitz's potential 
for simplicity. If we consider only the ground state interaction, last equation fits the exact potential 
within an error smaller than 0.1% for z > At- 

The ground and excited state potentials with the corresponding thermal contributions are 

V g (z, T) = V rr (z) + Vt (z) + V T (z) =: V (z, T) , (80) 
V e (z, T) = V rr (z) - Vt (z) - V T (z) = 2V rr (z) - V (z, T) , (81) 



^^For the Casimir-Polder limit, the condition Aq < z \t must be valid. 
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where V(z,T) is defined by equation (|80p . Considering the atom in thermal equilibrium with the 
field, we must consider the thermal average value of last potentials, since this average is the observed 
potential. If p is the probability of the atom to be found in its ground state, 1 — p is the probability 
of finding it in its excited state. In a thermal equilibrium at a temperature T, we have 



1 



P = 



1 _|_ c -fc Ar ' 

so that the average potential takes the form 

V (z, T) = tanh QkoAr) V (z, T) + J£ { *\ . (82) 

At low temperature, fco^T 1, one has 

V low (z, T) ~ V (z, T) - 2e~ koXT [v o fr (z) + V T (z) 

The average potential is given by the ground state potential minus the (fr) contribution weighted 
by the exponential factor e~ k ° XT -C 1, so that the average potential is practically due to the ground 
state potential only, as expected: the thermal photons do not have enough energy to excite the 
transition and populate the excited state, which makes negligible its contribution to the interaction. 
For example, at room temperature and for a transition frequency in the visible region, ko ~ 10/im -1 , 
the factor multiplying the (fr) term is about 7.1 x 10~ 34 . 

For the high temperature limit, however, things are diffent. In this limit, the excited state gives 
an important contribution to the interaction. In order to account the excited state influence on the 
interaction at high temperature, k Xx -C 1, it is necessary recalculate the integral in ([77)) with the 
help of Bernoulli's numbers, defined by the expansion [44 



1 2^ ^ 



e« - 1 2 s <H (2™) ! 

n— 1 v ' 

Using last expression in equation (|77|) . we obtain 



2n 



V T (z) = UA^ + hAP + TiAP , (83) 
hAft = / x 2 a' + (x) G (x) dx , (84) 



4lTZ 3 Jo 



hA^> = -—^ / x A a' + (x) G (x) dx = -Vf (z) , (85) 

1D7TZ 4 ' 







^ - Sll^ n j o ' X 2n+2 <( X )G( X )d X , (86) 

where 77 = Xt/2z and x = 2kz. Integrals in equation (|86p may be easily calculated with the help of 
equation (|56|) . which leads to 



hA^ = VJ r (z) 



coth [ —k^Xx 



k X T \ 2 



(87) 



The integral in equation ([M]) may be also calculated in a simple manner from equation (|56[) . but 
note that the function f(x) = x 2 G(x) is not analytical in x = and a direct application of (|56|) 
depends on the analyticity of /(x) over all complex plane. This problem may be avoided by splitting 
the integral in as follows 

f°° ( 1 n \ f°° sinx 

/ x 2 a' + (x) G (x) dx = nao ( xo sin xo — -x cos xq) —2 a' + (x) dx , 

Jo \ 2 J J a x 
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Figure 2: Dispersive potentials between a two- level atom and a perfectly conducting wall in a 
thermal bath at temperature T . We plot the temperature dependence of the thermal averaged 
potential V (z, T), as well as the potential associated to the ground state only Vlu (z, T), normalized 
by — Hu)oao/8z 3 (VJnean (0) and Vuf (#) respectively). The relative error |AV (0) \/Vlh {&) is also 
plotted. In the above figure it is assumed that z > At and 9 — 2ksT '/Tuoq. 



where xq = 2k$z. Last term in the previous equation can be calculated as 
f°° , sin a; f°° sin a; 1 f°° . , . 7rao 

/ a + \x) ax = oto I ax / a_ (x)sm.xax = — — (1 — cosa;o) ■ (89) 

Jo x Jo x x J 2 



Combining equations (f59"|) and equation ([54")) may be written as 



and the ground state potential ([50)) takes the form 

4 z 3 e fc A T _ ! ' 

Inserting last result into equation (182|) . the average potential may finally be written as 

V(z,6) = -^^f0tanh(l/0) , (92) 

where we defined the normalized temperature 9 = 2fc^T '/Tiljq — 2/ko\T- 

We derived equation (|92p by assuming the high temperature limit, which means 9 ~ 1 or higher. 
However, if we consider the opposite limit, 6 <C 1, we easily see that equation ([92)) reduces to 
equation (179)) , 

V (z, T) = (i _ 2e - fe ^ + 2e- 2fe °^ + ...) ~ _M!2° , 



15 



which is valid for all distances larger than the thermal length, that is, for z > At- Hence, equation 
(f9"2")l . which takes into account the excited state contribution, is also valid for all distances larger 
than the thermal length. 



E 

e 

E 

g 



Figure 3: A schematic spectrum of a system which simulates a two- level one and that can be used to 
check our result (|92p . Levels E e and E g are sufficiently close to each other so that thermal fluctuations 
can populate significantly both of them: E e — E g ~ ksT ~ 1/40 eV at room temperature. Level E c 
is so far from E e and E g (E c — E e 3> fceT) that, in practice, it can not be populated by thermal 
fluctuations. 

In Figure [5] we plot the potentials (|9"2")l and (|75|) normalized by the potential (f6"9"| . At low 

temperatures, as we have already anticipate, the excited state plays no role in the interaction and 
the ground state potential, which is given by Lifshitz's potential ([75]) . coincides with the average 
potential. However, at high temperatures, the difference between the Lifshitz's and the average 
potentials is evident, as can be checked in Figure[2]for values of 9 greater than 9 ~ 0.4 approximately, 
which is equivalent to a temperature about 2800-ftT x 1/im/Ao; at this point the relative error between 
the two potentials, defined as |AV|/VLif, is about 1.3% (for 9 = 1.0 the error is already ~ 24%). 

Nevertheless, the most interesting thing that happens as the temperature becomes higher is that 
the average potential saturate at a value independent of the temperature. In this case, the potential 
can be approximated by 

which coincides exactly with the potential (|69p in the limit of very high temperature, 9 — > oo. 

As mentioned before, expression (|69|) gives the vacuum term of the dispersive potential between 
the atom and the wall, for both ground and excited states of the atom, in the small distance limit 
(0 < z <C Ao). Further, when the system interacts with the field in a thermal state, thermal 
corrections to the potential become important only for distances larger than thermal length At- 
Equation (|93p is approximately valid for z > At at very high temperature, ^ > 1. This condition 
necessarily implies At <C Ao- Hence, combining last condition with the two just mentioned conditions 
for the separation z in which pontentials (|69|) and (|93p may be applied, one can see that London- van 
der Waals potential describes the interaction at very high temperature in the interval < z < oo. In 



Continuum 
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other words, for very high temperatures, the London- van der Waals potential is valid for all distance 
regimes. This is a quite curious result. 

However, though expressions (fuT))) and (|93[) are practically the same, there are some important 
differences in the physics behind them. The former is interpreted as the interaction between the 
instantaneous dipole and its mirror image so that retardation does not play any role at all on 
interaction and, as already mentioned, the main contribution is due to the radiation reaction. The 
second one is the net result of the competition of the fluctuations of the field acting on the ground 
and excited states of the atom simultaneously. 

From equations (|80ti81[) . one can see that thermal corrections have opposite signs and since at 
high temperature the probability of finding the system in its ground or excited states is practically 
the same, p ~ 1 — p ~ 0.5, the (/r) terms cancel each other and the temperature dependence tends 
to disappear in the average potential, as it occurs in (j9"3"| . We also see that in equations ([50)1 and 
(fSTj) the (rr) terms are the same, which leads to a net average potential different of zero. Hence, 
both equations (p)9"|) and (|93[) come from of the radiation reaction, though from different ways. The 
fact that retardation does not seem to have any effect on the average potential is a feature of thermal 
fluctuations which, for distances larger than At, destroy the influence of the finitcness of the velocity 
of light on the interaction, as one can see also in (l79|) . 

As a final comment, one may argue that equation ()93|) is an unrealistic one because there are 
not systems in the real world which can be treated as a two-level system at very high temperature. 
Actually, all known systems usually have an infinite number of energy levels and the coupling between 
them at high temperature should not be neglected, so that a different behavior of that preconized 
for the two-level system is expected [TB]. Eventhough, the two-level system model showed to be a 
very useful one in the study of the dispersive interactions and, particularly, in the comprehension of 
the roles played by (rr) and (fr) contributions. 

Though equation (|93[) is, in some sense, meaningless, the same can not be said to equation 
([92]) . It is plausible a system where the two lowest levels \g) and |e) are far enough to a third 
level \c), as sketched in the Figure [3l If the largest dimension Dq of this system is much smaller 
than the transition wavelength, Do <C litTicf (E e — Eg), and the temperature satisfies the condition 
E e — Eg ~ ksT -C E c — E e) then the equation (192|) will be valid and a deviation of the behavior 
predicted by equation (179")) should be observed for not too high values of temperature. 

4 Conclusions 

In this work we made use of the general expressions for level shifts obtained from master equation 
in order to study the dispersive potentials between a two-level atom and a perfectly conducting 
wall in the dipole approximation. We studied the potentials for the ground and excited states of 
the atom, which can be associated to the van der Waals and the resonant interactions, respectively. 
All distance regimes as well as the low and high temperature limits were treated. Considering 
the ground state potentials in the short distance limit, we reobtained the London-van der Waals 
potential given by (|69p , showing that it may be explained mainly by the reservoir reaction (radiation 
reaction) contribution to the level shift, as it is well accept [22]. For large distances, we rederived the 
Casimir-Polder interaction and showed that it can not be considered as a direct result of the vacuum 
fluctuations only, as usually intrepreted [42]. Otherwise, we showed that the field fluctuations and 
the radiation reaction contributions are of the same order and too much stronger than Casimir- 
Polder potential, but these effects "cancels each other" leaving a small observed interaction. For the 
case of the excited state potential, we emphasized that the maximum magnitudes of the potential 
occur approximately at the positions of resonant condition for stationary waves with wavelength Ao 
in the "cavity" defined by the atom and the wall. Taking into account the thermal corrections to 
the interaction, we rederived the short and large distance potentials at low temperatures, where the 
excited state contribution can be disregarded. However, for high temperatures, since the excited state 
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contribution should not be neglected anymore, we showed that its inclusion may cause considerable 
deviations from Lifshitz asymptotic potential. A curious result that we found is that at very high 
temperatures the potential saturate in a temperature independent value exactly equal to the London- 
van der Waals potential for all distance regimes. 

Acknowledgment 

The authors are grateful to Faperj and CNPq for the financial support. 



References 

[1] J.D. van der Waals, Over de continuiteit van den gas-en vloeistoftostand, Dissertagao, Leiden, (1873); 
J.D. van der Waals, The equation of state for gases and liquids, Nobel Lecture, (1910). 

H. Margenau and N.R. Kestner, Theory of Intermolecular Forces (Pergamon, New York, 1969). 

M. Reinganum, Ann. Phys., 10 (2), 334, (1903); 
M. Reinganum, Ann. Phys., 38 (8), 649, (1912). 

W.H. Keesom, Proc. R. Acad. Sci., Amsterdam, 23, 939, (1920). 

P. Debye, Phys. Zeits., 21, 178, (1920). 

R. Eisenschitz e F. London, Phys. Zeits., 60, 491 (1930); F. London, Phys. Zeits., 63, 245, (1930). 

Dieter Langbcin, Theory of Van der Waals Attraction, Springer Tracts in Modern Physics, Vol. 72, Springcr- 
Verlag, Berlin (1974). 

H.B.G. Casimir c D. Polder, Phys. Rev., 73 (4), 360, (1948); 

H. B.G. Casimir, Proc. K. Ned. Akad. Wet., 51, 793, (1948). 

J.E. Lennard-Jones, Trans. Faraday Soc, 28, 334, (1932). 
E. M. Lifshitz, Sov. Phys. JETP, 2, 73, (1956). 
E.S. Sabisky and C.H. Anderson, Phys. Rev. A, 7 (2), 790, (1973). 

I. E. Dzyaloshinskii, E.M. Lifshitz and L.P. Pitaevskii, Advan. Phys., 10 (38), 165, (1961). 

A.D. McLachlan, Proc. R. Soc. London A, 271 (134), 387, (1963); 
A.D. McLachlan, Proc. R. Soc. London A, 274 (AUG), 80, (1963). 

R. Kubo and K. Tomita, J. Phys. Soc. Japan, 9, 888, (1954); 
R. Kubo, J. Phys. Soc. Japan, 9, 570, (1957). 

G. Fcinbcrg and J. Suchcr, Phys. Rev. A, 2 (6), 2395, (1970). 

H. Morawitz, Phys. Rev., 187 (5), 1792, (1969). 

M. Babiker e G. Barton, Proc. R. Soc. London A, 326, 255, (1972). 

G. Barton, J. Phys. B, 7 (16), 2134, (1974). 

E.A. Power and T. Thirunamachandran, Proc. R. Soc. Lond. A, 401, 267, (1985). 
E.A. Power e T. Thirunamachandran, Phys. Rev. A, 50 (5), 3929, (1994). 

D. Meschede, W. Jhe e E.A. Hinds, Phys. Rev. A, 41, 1587, (1990). 

E. A. Hinds e V. Sandoghdar, Phys. Rev. A, 43, 398, (1991). 
W. Jhe, Phys. Rev. A, 43, 5795, (1991). 
W. Jhe, Phys. Rev. A, 44, 5932, (1991). 

H. Nha e W. Jhe, Phys. Rev. A, 54, 3505, (1996). 
G.H. Gocdccke and R.C. Wood, Phys. Rev. A, 60 (3), 2577, (1998). 
G. Barton, Phys. Rev. A, 64, 032102, (2001). 
G. Barton, Phys. Rev. A, 64, 032103, (2001). 

Y.J. Lin, I. Teper, C. Ching and V. Vulctic, Phys. Rev. Lett, 92 (5), 050404, (2004). 
M. Antezza, L.P. Pitaevskii and S. Stringari, Phys. Rev A, 70, 053619, (2004). 
E.V. Blagov, G.L. Klimchitskaya and V.M. Mostepanenko, Phys. Rev. B, 71, 253401 (2005). 

I. V. Bondarcv and Ph. Lambin, Phys. Rev. B, 72, 035451, (2005). 



18 



[33] G.L. Klimchitskaya, E.V. Blagov and V.M. Mostepanenko, arXiv: \quanTp h/0601032 vl, (2006). 

[34] T.N.C. Mendos and C. Farina, J. Phys. A, 39, 6533, (2006). 

[35] J. Dalibard, J. Dupon-Roc and C. Cohen-Tannoudji, J. Physique, 45, 637, (1984). 

[36] J. Dalibard, J. Dupon-Roc and C. Cohen-Tannoudji, J. Physique, 43, 1617, (1982). 

[37] C. Cohen-Tannoudji, J. Dupon-Roc and G. Grynberg, Atom-Photon Interactions: Basic Processes and Applica- 
tions, John Wiley and Sons, Inc., New York, (1992). 

[38] R. Kubo, Rep. Prog. Phys., 29, 255, (1966). 

[39] L.D. Landau e E.M. Lifshitz, Statistical Physics, Part I, Third Edition, Pergamon Press, Oxford, (1980). 

[40] G. Barton, Proc. R. Soc. Lond. A, 410, 141, (1987). 

[41] B.R. Holstein, Am. J. Phys., 69 (4), 441, (2001). 

[42] P.W. Milonni, The Quantum Vaccum: An introduction to Quantum Electrodynamics, Academic Press, New 
York, (1994). 

[43] National Institute of Standards and Tecnology, http://physics.nist.gov/PhysRefData/ASD/index.html (2006). 

[44] G. Arfken, Mathematical Methods for Physicists, 1st Edition, Academic Press Inc., New York, (1981). 



19 



